Effects of mixing in threshold models of social behavior 
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We consider the dynamics of an extension of the influential Granovetter model of social behav- 
ior, where individuals are affected by their personal preferences and observation of the neighbors' 
behavior. Individuals are arranged in a network (usually, the square lattice) and each has a state 
and a fixed threshold for behavior changes. We simulate the system asynchronously either by pick- 
ing a random individual and either update its state or exchange it with another randomly chosen 
individual (mixing). We describe the dynamics analytically in the fast-mixing limit by using the 
mean-field approximation and investigate it mainly numerically in case of a finite mixing. We show 
that the dynamics converge to a manifold in state space, which determines the possible equilibria, 
and show how to estimate the projection of manifold by using simulated trajectories, emitted from 
different initial points. 

We show that the effects of considering the network can be decomposed into finite-neighborhood 
effects, and finite-mixing-rate effects, which have qualitatively similar effects. Both of these effects 
increase the tendency of the system to move from a less-desired equilibrium to the "ground state" . 
Our findings can be used to probe shifts in behavioral norms and have implications for the role 
of information flow in determining when social norms that have become unpopular (such as foot 
binding or female genital cutting) persist or vanish. 
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I. INTRODUCTION 

In this paper, we investigate a simple model of behav- 
ior, the threshold model (TM) [U [2]. It consists of N 
individuals arranged in a network. Each individual Si 
{i = 1, . . . , N) has either adopted or rejected the behav- 
ior in question and has a tendency to switch to adopting 
(rejecting) if the proportion of individuals in its neigh- 
borhood adopting the behavior is greater (less) than its 
(constant) threshold T^. Individuals are chosen at ran- 
dom to be "updated" - i.e., to consider, and possibly 
change ("flip") their state. We make an analogy with 
physics by thinking of the individual's state as a 'spin' 
with value +1 (—1) for those who adopt (reject) the be- 
havior. 

Threshold models are relevant to questions of how pat- 
terns of behavior persist, even when attitudes change, 
and how these patterns can sometimes change rapidly. 
A currently relevant example is the practice of female 
genital cutting (FGC), which goes back at least to an- 
cient Egypt [3^. Despite a public health consensus that 
the practice is harmful [4], traditional practice remains 
widespread in various societies [SJ [5]. A similar exam- 
ple is the Chinese practice of footbinding, which was 
widely practiced for hundreds of years, before disappear- 
ing rapidly ^ . These practices can be considered in the 
context of the theory of "social norms" , behaviors which 
individuals prefer to follow, given that they think that oth- 
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ers will conform,, and that others expect them to conform 

Many similar individual-based models are also based 
on individuals making binary choices [H [10] . Usually, 
the voter model [TT] is associated with imitation process, 
since a randomly chosen individual adopts the behavior 
of one of its neighbors. In this sense, the TM puts the so- 
cial pressure in the framework: the adoption or rejection 
of the behavior by an individual depends on the current 
level of adoption in its neighborhood [T5]. The majority 
rule model (see [13]) is a special case of the TM, since all 
thresholds are one half (the randomly chosen individual 
tends to flip if the majority of its neighbors have opposite 
spin). 

It has been shown that the majority rule model can 
be described by the classical Ising model with zero ex- 
ternal magnetic fleld |T3j and that the general TM can 
be described as a random-fleld Ising model (RFIM) |16) . 
The study of RFIMs in physics often focuses on critical 
temperature phenomena |17| or metastable states and 
hysteresis loop phenomena at zero temperature [THl [IH] ■ 
Instead of using the notion of the thermodynamic tem- 
perature, where individuals probabilistically flip in a non- 
preferred direction, see, for example, [201111], we chose 
to set the thermodynamic temperature to zero and study 
the effects of mixing on the dynamics. 

We simulate our model on a two-dimensional lattice, 
with global mixing. We implement mixing by allowing in- 
dividuals to exchange places within the network at rate 
/i (relative to the update rate) . The importance of mix- 
ing in sociological and ecological studies has been demon- 
strated in other contexts [^^H^ . Introducing global mix- 
ing on a two-dimensional lattice is similar in concept to 
using a "small- world" network [25] . Both cases have reg- 
ular connections, and random global connections - the 



difference is that we implement random global connec- 
tions by switching individuals. 

We simulate behavior change by either choosing an in- 
dividual at random to update or mixing two individuals 
in each step of the simulation. Mixing consists of ex- 
changing two randomly chosen individuals rather than 
updating in a given simulation step, with probability 
fi/{l+fi), so that we have an average of /i switches per up- 
date. We increment the clock by 1/A^ per update event. 
This gives us update events at rate 1 per individual and 
mixing events at rate /i per individual. When we mix in- 
dividuals, we exchange their states and thresholds, leav- 
ing the network otherwise unchanged. A synchronous 
process or a pure Poisson process would be expected to 
give qualitatively similar results, but this asynchronous 
process is simpler to simulate, and can be analyzed us- 
ing an ordinary differential equation (ODE) framework 
derived from master equations. 

Most analytical results in the field of TMs/RFIMs are 
obtained using mean-field approximations [271 128] . This 
can be achieved either by considering a complete network 
(where every node is a neighbor of every other node) , or 
by setting the mixing rate /i ^ 1. 

The intermediate mixing case /i ^ 1 is not so easily 
treated. If we write equations for the moments of differ- 
ent order for the distribution of states among individuals, 
we get a hierarchical system of coupled equations p8] . 
There are then various methods to "close" the system 
by approximating higher moments in terms of lower mo- 
ments [29H31j . 

In [32], the authors considered a MR model, and con- 
cluded that the behavior of their system resembles the 
movement of a Brownian particle in a potential field that 
is unknown a priori. We describe such an "effective po- 
tential" function for our threshold model and calculate 
the analytical potential form for the mean-field version 
of the model. 

We can use the effective potential to provide an ad- 
ditional perspective on the dynamical properties of the 
system. The bifurcation where the system changes from 
having one stable equilibrium to two, for example, cor- 
responds to a change from a single-welled to a double- 
welled effective potential function. In terms of an Ising- 
like model, this would correspond to a phase transition 
of the first order. This bifurcation is relevant from a so- 
ciological point of view, since a transformation from a 
potential consisting of two wells to a potential consisting 
of one well, due to a change of mixing rate, could give 
rise to sudden abandonment or adoption of a social norm. 
In contrast, if such transformation does not occur, even 
when one well is much deeper than the other, this might 
help to explain why a human society sometimes contin- 
ues to support a fairly unpopular social norm for many 
years [8]. 

In this paper we explore the dynamics of this system 
using a Gaussian threshold distribution. We expect other 
flexible distributions to give similar qualitative results. 
For example, even uniform distributions show the same 



basic bifurcations that we explore here \T. We have also 
tried simulations with bimodal superpositions of Gaus- 
sians, and again seen qualitatively similar results. 

Here is how we organize the rest of the paper. First we 
consider the mean-field dynamics which are given by the 
fast-mixing and large-scale {N S> 1) limits. Then we con- 
sider the intermediate mixing rates and state the main 
results of our paper: we discuss the bifurcation phenom- 
ena found in the TM and we demonstrate the appear- 
ance of a manifold in the dynamics that is approached 
by any trajectory of the TM. Finally, we interpret our 
main results from a sociological point of view, and draw 
conclusions. 



II. MEAN-FIELD APPROXIMATION 

First, consider the case in which the neighborhood size 
is so large that each spin is connected with all other spins 
in the network. In this case, we recover Granovetter's 
threshold model for collective behavior [2] , with dynamics 
in which the probability of a spin being updated from 
minus to plus is given by P(t \y) = (1 — y)F{y), where 
y denotes the proportion of plus spins and F{-) is the 
cumulative distribution function of the thresholds' PDF 



F{x) 



/(OdC 



(1) 



The probability of a spin being fiipped in opposite di- 
rection from plus to minus is P(| \y) = y{l — F{y))- In 
this case, master equations can be written in terms of the 
probability function p{yk,t) {yk ~ k/N), which provides 
the probability to find the TM in a state with k spins in 
a plus state and N — k spins in a minus state at a given 
moment of time t: 



dp{yk,t) 



dt 



= m\yk-i)p{yk-i,t) 



-V{i\yk+i)p{yk+i,t) 



nt\yk)p{yk,t), 

k = 0,...,N, 



(2) 



where P(t \yk) = P(t IVk) + P(i \yk). 

Letting N ^ oo and scaling the time as i — >■ Nt, we 
can treat the discrete variable yk as continuously chang- 
ing y g [0, 1] and transform (pi) to the Hamilton-Jacobi 
equation, which is the first order partial differential equa- 
tion: 



dpjy, t) 
dt 



d_ 
dy 



[{F{y)-y)p{y.t)] 



(3) 



During such transformation, the diffusive terms, consist- 
ing of second order partial derivatives, vanish due to the 
large-scale limit TV ^ 1. 

If the initial configuration is strictly defined such that 
p{y,0) — 5{y — yo), where S{-) is the Dirac delta, the 




FIG. 1. Simulations of the threshold model (TM) on a two- 
dimensional lattice of size 100^ with no mixing among individ- 
uals, and different neighbourhood sizes. The thresholds are 
normally distributed with the mean 0.45 and standard devi- 
ation 0.3. The initial pattern of thresholds and initial states 
is the same for all simulations shown. Activated individuals 
are shown in color 



solution of (Is]) is given by the following ODE, see p. 53- 
54 M: 



y = ^ = F{y) " y, y(o) = 2/0 



(4) 



The equilibria of this system are all values y* for which 

P{y*) = y*- 

Notice that Q can also be written in terms of the 
potential function: V{x) = — J (F(^) — ^)d^, such that 
-^ = —V'{y). Thus, the equilibrium points can be also 
defined by the extrema of V{y). 

We now consider a case where each individual's up- 
dates depend on the states in a finite neighborhood. First 
simulations of the TM on a two-dimensional lattice with 8 
nearest neighbors for each individual and with no mixing 
among them reveal complex patterns, shown in Fig. IT] 
This figure presents initial, intermediate and final states 
of the lattice for four different neighborhood sizes, but 
with identical initial distributions of states and thresh- 
olds. We see that increasing neighborhood size can shift 
the outcome of the system's dynamics from a low level 
of conformity to a very high level. Moreover, the equi- 
librium distribution preserves some noticeable clustering 
for small neighborhood sizes. 

However, in the large-scale {N — > oo), fast- mixing 
(/i — ^ oo) limit, the behavior of the TM can still be 
described analytically by a mean-field model, in which 
a spin and all its neighbors are chosen de novo at each 
update event. 

The probability that a randomly selected individual 
will choose to adopt is equal to the probability that the 
activation level of its randomly selected neighborhood 






FIG. 2. (Color online) Functions Fn{y) and corresponding 
potentials Vn ijj) (in the inset) with respect to different neigh- 
borhood sizes: n = A (magenta), 12 (red) and 24 (green). The 
dashed blue curve corresponds to the case of infinitely large 
neighborhood size, Foo{y) = F{y). The thresholds' PDF is 
Gaussian with the mean 0.6 and standard deviation 0.225. 
The crosses show the results of numerical simulations of the 
TM on a two-dimensional lattice of the size 100^ at mixing 
rate jj, — 4 and j/(0) — 1 

exceeds its threshold. In a regular network where each 
individual has n neighbors, this is given by: 



fc=0 



F„iy)^Y.F[^)^ny'{i-y) 



i—k 



(5) 



where C,^ — j^u^Lkv is a binomial coefficient (see [T6]). 
Then (|2])-(|4]) remain valid for this system once we substi- 
tute F^y) for F{y) in ^, and they give us the dynamics 
of the TM with finite neighborhood size in the mean-field 
approximation. (This approach also works for networks 
of variable degree; if neighborhood sizes are distributed 
with probability density 7'(n), we average over the dis- 
tribution to get F-p = i2n=o'^i^)Fniy)-) 

Fig. [2] illustrates the difference between the functions 
([5]) and (IT]) as well as the difference in the mean-field 
dynamics of the TM for different neighborhoods. Pa- 
rameters of the thresholds' PDF / for Fig. J2] were chosen 
in such a way that there are two stable equilibria for 
large neighborhoods, and only one stable equilibrium for 
small neighborhoods. Different neighborhood sizes can 
lead to very different outcomes, even when distributions 
and initial conditions are the same. For example, in sim- 
ulations starting with everybody adopting the behavior 
(j/(0) = 1), the TM reaches a high equilibrium (few indi- 
viduals change), when neighborhood size is large, but a 
low equilibrium (almost everybody rejects the behavior) 
when neighborhood size is small. Note that simulations 
done on a two-dimensional lattice of size 100^ at mixing 
rate /i = 4 give a good approximation to the large-scale, 
fast-mixing limit in this case; later we will show that this 
is not true for smaller mixing rates, though. 

In case of a Gaussian distribution for the thresholds, 
the curve y' = Fn{y) has up to three crossings with 
the diagonal y' = y. If there is only one crossing with 
the diagonal, there exists a globally stable equilibrium 
y- e [0,1]. If there are three crossings, we have three 



equilibrium points, which we denote as y-, y* and y+, 
such that y- < y* < y+- Two of them, y±, are stable 
equilibria and one of them, y^ , is an unstable equilibrium. 

We can define a potential, analogous to the mean-field 
case: Vniy) = /^(^n(0 ~ C) ^C- Then y± are the min- 
ima and y* is the maximum of Vn{y), see Fig. ^ (inset). 
The case of two crossings represents the bifurcation point 
between these two generic cases. 

If the mean of the Gaussian distribution is not exactly 
0.5, the potential function Vn{y) is asymmetric, with one 
well deeper than the other. Without loss of generality, 
we assume that the norm is intrinsically unpopular (i.e., 
the mean of the threshold distribution > 0.5), so that 
the "lower" equilibrium y_ corresponds to the deeper 
well, and the "upper" equilibrium and ?/+ to the shallower 
well, when it exists. These values refer to the case where 
/i -> oo. For clarity, we will sometimes add oo as a 
superscript. 



III. INTERMEDIATE MIXING 

Simulations show that reducing the mixing rate away 
from the fast-mixing limit has a similar qualitative effect 
to reducing neighborhood size (as seen in Fig. [2]) . In the 
case where the mean-field system has one stable equilib- 
rium, reducing mixing rates does not lead to a qualitative 
change in the dynamics. In the case where the mean-field 
system has two stable equilibria, as the mixing rate gets 
smaller we often find a bifurcation to a single equilib- 
rium; i.e., the equilibrium with the shallower potential 
well disappears. 

We consider a two-dimensional lattice, with initial ac- 
tivation level y(0). If we simulate, starting from a value 
between the two stable equilibria, the system will move 
to the upper equilibrium with probability p_(_; otherwise 
it moves to the lower equilibrium. The result depends on 
the random selection of thresholds, initial states and the 
order in which sites are updated. 

Fig. [3] shows how the probability p+ depends on the 
mixing rate /i, using two different scaling approaches. In 
either case, as we move away from the fast-mixing case 
(from right to left on the figure), the system becomes 
increasingly certain to end in the deeper well, and even- 
tually the shallower well disappears altogether. 

The main picture of Fig.lSlshows the probability p+ vs. 
IJ, for values of y(0) > y^ . In this case the system stops 
in the shallow well for large values of /x, and moves to the 
deeper well for smaller values. This transition becomes 
steeper as the size of the network increases. The curves 
for a given starting point intersect where p+ = 1/2. That 
is to say, for a given value of fi, the value of y(0) that falls 
"in the middle" of the two wells - so that the system is 
equally likely to go to either one - does not change with 
lattice size. We call this value y'^, because it is related 
to the unstable equilibrium y^, but not equivalent (as 
we will see below). The dependence of y^ on fi has a 
hyperbolic shape and its minimal value p, is reached at 
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FIG. 3. (Color online) The probability p+ that the TM will 
approach the upper equilibrium as a function of /i or scaled 
fj,/L (the inset). The TM is on a 2d-lattice of size N = L^ 
with the neighborhood 8. The thresholds' PDF is Gaussian 
with the mean 0.55 and standard deviation 0.225. Different 
colors stand for different initial values j/(0), while the symbols 
stand for different values of L (see the legend in the bottom 
right corner). The bifurcation value jl{y{Q) = 0.69) « 1.556 
is indicated. To estimate the probability, we performed 10"" 
simulations with different random initial individual states and 
thresholds, and update order. 




FIG. 4. The dependence of j/^ on A*- The TM is posed on a 
2-d lattice of size 200^ with neighborhood 8 and Gaussian dis- 
tribution of thresholds with the mean 0.55 and standard de- 
viation 0.225. The extrapolation curve, shown by the dashed 
line, gives Ji — niy'^ = 1.0) ^ 0.168. 



j/x = Ij which is shown in Fig. Hi 

The inset in Fig. [3] shows the same data, with p+ plot- 
ted against a scaled version of the mixing rate n/L (where 
L = N^/'^ is the length of the two-dimensional lattice). 
A surprising pattern emerges. For y(0) = yj°, all of 
the curves approximately align onto a single curve, with 
p+ —^1/2 for ^ — )> 00, as we expect, since we are ap- 
proaching the well-mixed case, where y(Q) is the unstable 
equilibrium. For other values of y(0), the curves do not 
intersect in this scaling: instead, as N gets larger, the 
system becomes less likely to "switch" to the equilibrium 
on the other side of y"^ . 

If we visualize the trajectories in the phase subspace 




FIG. 5. The projection curve F^ on a two-dimensional lat- 
tice of size 800^ with neighborhood 8 and the Gaussian dis- 
tribution of thresholds (the mean 0.55 and standard devia- 
tion 0.225 and mixing rate equals 0.278. The trajectories, 
shown in red, start from y(0) = i/10 (i — 1, ... ,9), while the 
green curve consists of the trajectories initiated at y{0) = 0, 
y(0) = 1 and t/(0) = Vx = 0.76. Initially, the states and 
thresholds are distributed randomly among individuals, hence 
all initial points fall along the curve F„{y) = F^(y) (dashed 
curve) 
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FIG. 6. (Color online) The mean transition time (Tr), neces- 
sary for the TM to evolve from complete activation y{0) = 1.0 
to a given intermediate value 2/(Th) = yT ~ 0.6752. The TM 
is posed on a two-dimensional lattice of size N — L^ with 
neighborhood 8 and Gaussian distribution of the thresholds 
with the mean 0.55 and standard deviation 0.225. The inset 
illustrates the distribution of Tn at jj, — 0.1 for L = 100 (ma- 
genta), 200 (green) and 400 (blue) 



{y,dy/dt), we find that all of them approach the same 
curve F^, shown in Fig. [51 presumably because they are 
collapsing onto a lower-dimensional slow manifold. Thus, 
the behavior of the TM can be well-approximated by the 
ODE: dy/dt — F,l^{y) — y, on some time interval t e 
[ii,oo), which is similar to (B, where F^ = F„. The 
equilibrium points yt^ can be determined as F^(]jl^) — yt 
and the effective potential can be introduced by V^{y) = 
~ I^i^niO ~ C)dC- Thus, we can describe the behavior 
of the TM qualitatively, by studying the properties of the 
manifold-projection curve Fj^. 

When y approaches or 1, mixing does not affect the 
dynamics. Therefore, we can construct at least part of 
the curve F^ (for any value of /i) by simulating trajec- 
tories starting from y{0) — 1 and y{0) = 0. When there 
is only one equilibrium in [0,1], this method generates 
the whole projection curve. When there are two sta- 
ble equilibria, this method generates only the part "out- 
side" them. Completing the curve requires that we start 
simulations from one or more intermediate initial points 
y{0) e (0, 1). In fact, we need only one additional start- 
ing point, which is precisely y^, since any trajectory ini- 
tiated at that point goes through the point yf^ on the 
curve F^ to the upper or lower equilibrium with equal 
probability one half. 

Note that if we take the minimal value ft of the mix- 
ing rate, which can defined using Fig. l4J F^ will have 
two fixed points, one of them will be a double root of 
^niy*) = y*^ which corresponds to the bifurcation, de- 
scribed above. 



Transition times 



For mixing rates /z < /2, the system will always move 
towards the lower equilibrium y'^. We therefore explore 
the "transition time" Tr - how long it takes to move from 
the fully activated state y{0) = 1.0 to some intermediate 
activation level, chosen to be near, but to the right of, 
the lower equilibrium. 

From simulations, we see that the distribution of tran- 
sition times becomes narrower for larger N, see the in- 
set Fig. [6] Hence, it is important to look how the mean 
value (Tr) changes for TV ^ 1. It turns out that the 
dependence (Tr) vs. fj, is concave and has a minimum 
at some intermediate value ji. We see that it becomes 
arbitrarily large for small mixing rates and exponentially 
decreases as n approaches jl. From the other side the 
value (Tr) increases as /i becomes closer to ft, see Fig. |6] 

The minimum in the transition-time curve can be ex- 
plained in terms of two countervailing effects of increased 
mixing. When the mixing rate is very slow, any changes 
in behavior take a long time to spread through the lat- 
tice. When the mixing rate is high, for our parameters, 
exchange between neighborhoods tends to preserve the 
"upper equilibrium" , leading to an exponentially slow 
transition on the finite lattice (and no transition at all 
for the infinite system). Thus, the most rapid transition 
from the activated state to the lower equilibrium occurs 
at an intermediate value. 



IV. CONCLUSIONS 

Offering an explanation why collective behavior can 
shift abruptly from avoidance to adoption of an alter- 
native, or vice versa, [2] provides an example where 
a slight change in distribution of individual thresholds 
leads to completely different outcomes: a system initially 
at one stable equilibrium switches to another one due to 
a change of the threshold of one individual. This change 
can be visualized as a change in the shape of the -F„- 
curve, which we will call the "activation curve" . We in- 
vestigate other factors that can change this curve and 
lead to similar phenomena, including abrupt changes in 
outcome when an equilibrium disappears. 

We model a population on a lattice, with a finite inter- 
action neighborhood, and "mixing" - implemented by ex- 
changing random individuals. To disentangle the effects 
of neighborhood size and "locality" , we first considered 
a lattice with finite neighborhoods in the infinite-mixing 
limit. We showed that the effect of finite neighborhoods 
is to "flatten" the activation curve, often leading even- 
tually to the elimination of the "weaker" equilibrium, as 
neighborhood size gets smaller. 

We then consider the effect of "locality", by reducing 
the mixing rate to be of the same order as the update 
rate. This system is harder to analyze, but we show that 
it tends to converge towards a manifold, whose projec- 
tion can be interpreted in a way similar to the activation 
curve. This interpretation allows us to define an effec- 
tive potential in the finite mixing case, which can aid in 
qualitative analysis. We find that the effect of locality 
on the projection curve is similar to the effect of finite 
neighborhood size: it flattens the curve and eventually 
leads to the disappearance of the weaker equilibrium. 

The flattening due to finite neighborhood size Fig. [2] 
can be understood in terms of averaging. If each indi- 
vidual evaluates a random, finite subset of the popula- 
tion when updating, the realized activation curve is a 
weighted average of the original curve. This averaging 
tends to flatten out curvature: in the limit of considering 
a single neighbor, the activation curve becomes a straight 
line. Finite mixing has a similar effect Fig. [5] Individu- 
als' states will be correlated with those of their neighbors, 
since they are responding to each other. This increases 
the variance in neighborhood activation perceived by in- 
dividuals, for a given value of the mean activation, ac- 
centuating the effect of averaging and further smoothing 
the activation curve. 

Here we, as others in the past |TnillH]i use the mapping 
between the threshold model and the random field Ising 
model so that it is possible also to apply tools from sta- 
tistical physics to the question. Another possibly useful 
analogy can be made between the TM and a spin gas. 
We neglect the structure of the network and consider 
particles stochastically moving in uniform medium, and 
affected primarily by nearby particles. In this case, the 
mixing rate can be associated directly with thermody- 
namic temperature. There is then an analogy between 



the tendency of all spins to be at the lower equilib- 
rium for small mixing rate in the original system, and 
low-temperature Bose-Einstein condensation in the spin 
gas [33] . This mapping may be worth future study. 

From sociological point of view, mixing is associated 
with the rate of information flow in a given society or 
people mobility. We might imagine "activists" who have 
high mixing rates, and who are eager to change the preva- 
lent behavior. We have seen that sufficiently large mix- 
ing rates actually prevent the system from switching to a 
desirable equilibrium, so that an unpopular social norm 
persists, while low mixing rates facilitate the abandon- 
ment of the social norm. However, we should be aware 
that as mixing rates are reduced the transition can be- 
come very slow, so that a situation in the middle, when 
the mixing rates are not large and not small enough, may 
be the most preferable. 
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APPENDIX A: ISING MODEL FRAMEWORK 

In our model, each node i G {1 . . . N} of the network 
has a state (or spin) Si and a (constant) threshold Tj. 
The classical Ising model translates to the a majority 
rule (MR) model, where all thresholds are exactly 0.5: 
a spin tends to flip to the state where it will be aligned 
with more than one half its neighboring spins, see Ch. 8 
|28| . This system is Hamiltonian with the energy func- 
tion: n = -hJ2^■,J 

refers to the network neighbors of node i. 

When the thresholds are randomly distributed with a 
given probability distribution function (PDF) , the model 
becomes equivalent to the spin system under a magnetic 
field which describes effects of locality between spins, and 
such that the strength of nearest interactions depends on 
the connectivity of the network. In this case, the system 
also obeys Hamiltonian dynamics and its energy function 
has the form 



^i;jg(j) SiSj, where i € 1 . . . iV, and (i) 



n 



E 



SiSj 

2ni 



J2{2T,-l)s^ 



where n^ is the number of connections for a spin Si and 
Ti is a given threshold of it. Thus, the induced magnetic 
field is hi=2Ti-l. 

To simulate the TM dynamics, the following underly- 



ing update rule is posed for each update event 



Si i-> sign - 



dn 



= sign 






Sj - h. 



(6) 



while the thermodynamic temperature, determining the 
rate of random flips of spins, is set to zero. Hence, we use 
only the first part of the Metropolis algorithm [3S] that 
consists only of ^ in order to simulate the dynamics of 
the TM. The second part when the spin might be flipped 
even if it was not updated due to (|6| is omitted. 

Note that ^ indeed allows a sociological interpreta- 
tion of the TM: if the proportion of plus spins (individ- 
uals adopting the social norm) in the neighborhood of a 
spin Si is written as 1/° = ^'^j^u) 2^ , such that rii 
is the connectivity of a spin Si, then (|6| transforms to 
the following form: Si 1— > sign(2/° — Ti). In a particu- 
lar case of the MR, it translates to the simplest form: 
Si ^ signJ2je{i} «J = sign(2/° - 1/2). 



APPENDIX B: TECHNICAL DETAILS OF 
SIMULATIONS 

In our simulations, individuals' initial states and 
thresholds are independently identically distributed with 



a given initial activation level y(0) and PDF of thresh- 
olds. After that, we initialize the simulation process, us- 
ing the Metropolis algorithm. At each event, we either 
update with probability (1 -I- /i)~^ a state of a randomly 
chosen individual or exchange with complementary prob- 
ability the locations of two randomly chosen individuals. 
We associate the time only with update events, by defin- 
ing the time quanta 1/N. We say that the TM is located 
near the equilibrium point if it fluctuates near it over a 
sufficiently long period of time, compared with the time 
of observation. 



To construct Fig. [3] and Fig. [4J we considered the TM 
with i/f° < y(0) < 1/°^ . In this case, we expect the system 
to move toward either the lower or upper equilibrium. 
We run the simulations until they traverse 85% of the 
distance from the starting point to one of the mean-field 
equilibria; p_|_ is the probability that it has moved toward 
the upper equilibrium. 
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